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Abstract 

The influence of a single magnetic cluster in a non- magnetic host metal on the spin current j^'*^ 
and the charge current j in the vicinity of a ferromagnetic STM tip is studied theoretically. Spin- 
flip processes due to electron interaction with the cluster are taken into account. We show that 
quantum interference between the partial waves injected from the STM tip and those scattered by 
the cluster results in the appearance of components perpendicular to the initial polarization of the 
spin current which obtain a strongly inhomogeneous spatial distribution. This interference 
produces oscillations of the conductance as a function of the distance between the contact and 
the cluster center. The oscillation amplitude depends on the current polarization. We predict a 
strong magneto-orientational effect: the conductance oscillations may grow, shrink, or even vanish 
for rotation of the cluster magnetic moment by an external magnetic fleld. These results can 
be used for the determination of the /Xgg for magnetic clusters below a metal surface. 

PACS numbers: 61.72Ji, 73.40Cg, 73.63Rt, 74.50+r 
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I. INTRODUCTION 



Subsurface defects, such as impurities and vacancies, result in oscillations of the conduc- 
tance as a function of the position of a Scanning Tunnelling Microscope (STM) tip relative to 
the defect position (see, for example, [H, 0, 0]). These Friedel-like oscillations originate from 
interference between electron waves directly transmitted through the contact and waves that 
are scattered by the defect and reflected back by the contact. The theory of STM conduc- 
tance in the presence of a single non-magnetic point-like defect below a metal surface in the 
vicinity of the contact has been developed in Refs. [1, Hj. First the signature of Fermi-surface 
anisotropy in STM conductance in the presence of defects had been analyzed theoretically 
in detail in Ref. [gI]. In the paper [3| the general results of Ref. [i] was applied for the the- 
oretical investigation of the conductance of a tunnel point contact of noble metals in the 
presence of a single defect below surface. A pattern of the conductance oscillations, which 
can be observed by the method of scanning tunneling microscopy, was obtained for different 
orientations of the surface for the noble metals 0]. Recently it had been confirmed experi- 
mentally that Fermi surfaces can be imaged in real space with a low-temperature scanning 
tunneling microscope when subsurface point scatterers are present [sf . The effect of Kondo 
scattering by a subsurface magnetic impurity on the conductance of a tunnel point contact 
has been analyzed theoretically in Ref. j^. 

The applicability of STM can be extended for the study of magnetic objects below the 
surface of a conductor when a magnetic material is used for the STM tip such that the 
electric current is spin polarized (SP) (for review of SP-STM see [13] )• An important feature 
of SP-STM is that the spin polarized current influences a magnetic object in a non-magnetic 
matrix, producing so-called spin transfer torque (for review, see For example, near 

a point contact, where the current density is large, the spin transfer torque can be strong 
enough to reorient the magnetization of ferromagnetic layers in magnetic multilayers (l2| . 
Such investigations are very important for the development of innovative high-density data 
storage technologies. 

In this paper we consider theoretically the conductance of a tunnel point-contact between 
magnetic and non-magnetic metals in a SP-STM geometry. A magnetic cluster is embedded 
in the non-magnetic metal in the vicinity of the contact. The changes in the spin-polarized 
current and the spin transfer torque that influences the magnetic moment of the cluster are 
analyzed. We study the dependence of the amplitude of the conductance oscillations as a 
function of the STM tip position, on the relative orientation between the spin polarization 
of the tunnelling electrons and the magnetic moment of the cluster ^tg^. 



II. MODEL OF THE CONTACT AND BASIC EQUATIONS 

Let us consider a tunnel contact between a semi-infinite half-space 2; > of a nonmagnetic 
metal (the sample) and a sharp tip of a magnetic conductor (Fig. 1). A voltage V is applied 
between the tip and the sample. The electrical (and spin) current flows through a small 
region of the surface at z = near the tip apex that is closest to the sample. This system 
models the geometry of a SP-STM experiment. The tip magnetization (in real SP-STM 
the magnetization of the last atom l3]), which we choose along the z-direction, deflnes the 



direction of the polarization of the tunnel current. Such magnetization can be obtained 



for example, for a Fe/Gd-coated W STM tip [13|. In the non-magnetic metal we place a 



spherical single-domain magnetic cluster having a magnetic moment (Fig. 1). As flrst 
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FIG. 1: Model of the contact in a SP-STM configuration, with a subsurface magnetic cluster near 
the tunnel contact point. Spin-polarized electrons tunnel into the non-magnetic metal in the small 
area below the STM tip. The trajectories of electrons that are scattered by the spherical magnetic 
cluster are shown schematically. 



predicted by Frenkel and Dorfman 13] particles of a ferromagnetic material are expected to 
organize in a single magnetic domain below a critical particle size (a typical value for this 
critical size for Co is about 35nm). Depending on the size and the material, the magnetic 



moments of such particles can be /igg ~ 10 — 10 /ig, where /ig is the Bohr magneton [15 



Below, we only consider electron scattering by the magnetic cluster, assuming the mean free 
paths for all other processes (electron spin- diffusion length, electron-phonon mean free path 
and others) are much larger than the distance between the contact and the cluster center 

In order to describe the spin-polarized electron states of our system we use the approach 
proposed in the works of Slonczewski and Berger [l^, 21 1: all calculations are performed 



by means of independent single particle spinor wave functions \l/(ri; cr^) of electrons with 
opposite spin directions, where r/ and ai are the position vector and the spin direction, 
respectively, for each spin orientation / = 1,2 . We use the representation ai^2 =T, i, in 
which the spin projection on the z axis, Sz = ±|, is used. This approach corresponds 
to reducing the many-particle problem of a partially polarized electron system with non- 
zero average spin to a double-particle problem for electrons in a pure (completely polarized) 
spin state. Neglect of electron-electron interactions enables us to separate the double particle 
Schrodinger equation into two independent equations for \l/(r;; o";). In our case this separation 
is valid if fi^g ^ /Xg, i.e. the electron-electron exchange interaction is negligible compared 
to electron exchange interaction with the cluster. Generally, the moment of the cluster 
in a non-magnetic metal takes an arbitrary direction. This direction (the angle a in Fig. 1) 
can be held fixed by an external magnetic field H^, the value of which is estimated as 
Ha — T/fi^ff, where T is the temperature (see, for example, Ref. [l3])- For /Xgg ~ lO^/Xg 
and T ~ IK the field is of the order of O.IT. We assume that is much smaller than 
the magnetocrystalline anisotropy field of the STM tip, i.e. the direction of the external 
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magnetic field controls the direction of the cluster magnetic moment but its influence on the 
spin-polarization of the STM current is negligible. If the magnetic moment f-i^g of the cluster 
is 'frozen' by the field Ha the problem becomes a stationary one. Also we take the distance 
between the contact and the cluster tq to be much smaller than the radius = cpp/eHa of 
the electron trajectory in the applied external magnetic field H^, and we do not take into 
account trajectory magnetic effects, which have been analyzed in Ref. [l8|. The condition 
'^o th together with inequality fi^Ha/e <^ 1 (e is the electron energy) allow neglecting 
the magnetic field in the single-electron Hamiltonian and considering as an independent 
external parameter. 

The external magnetic field may result in a modulation of the tunnel current due to 



electron spin precession [19|. For our problem such precession would become noticeable 
when the transit time for the electron motion from the contact to the cluster t ~ vq/vf 
is larger than 1/ujl, where vp and are the Fermi velocity and Larmor frequency. The 
inequality mentioned above, -C th, is equivalent to the condition ui^t <C 1, so that the 
effect of electron spin precession is negligible. 

Thus, under the assumptions outlined above the spinor wave functions \I'(r/;cr/) satisfy 
the one-electron Schr odinger equation 

Vf - s) m{ri; ai) = -f/(rO$(r,; a,), (1) 



2m* 



where m* is the effective mass of the electrons, and U{r) is the interaction potential of the 
electrons with the cluster. The matrix U consists of two parts, describing the potential as 
well as the magnetic scattering, 

Uir) = (gT+ -Lj ^^^a] Di\r - ro|), (2) 



where g is the constant describing the non-magnetic part of the interaction (for g > the 
potential is repulsive), J is the constant of exchange interaction, /x^g = ^^^(sina, 0, cosa) 
is the magnetic moment of the cluster, a = (cxjO'yjO'z) with the Pauli matrixes, and / 
is the unit matrix. Z)(|r — ro|) is a spherically symmetric function localized within a region 
of characteristic radius ro centered at the point r = Tq, which satisfies the normalization 
condition 

^c/rD(r') = l. (3) 



Equation ([T]) must be supplemented with the common boundary conditions for continuity 
of the wave function and its normal derivative on the metal surface. 

We assume that the potential U (r) is small and use perturbation theory. As a first step 
we should find the solutions \l/^°)(r;; cr/) of Eq. ([1]) for U{r) = 0. Generally, this solution 
depends on the model chosen to represent the tunnel barrier. For any suitable model for 
the potential barrier the wave functions '^^'^\ri; ai) describe the spreading of the electron 
waves in the metal from the small contact region on the surface. Here, we use the results 
of Refs. 20|, in which the tunnel contact is modelled in the form of a circular orifice of 



radius a, with a large amplitude potential barrier UqS{z). 

In order to describe the spin polarization of the STM current we introduce different 
magnitudes for the wave vector k^- for spin-up and spin-down electrons (for the same energy 
e) before tunnelling from the tip. This difference results in different amplitudes t(ko-) of 
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the electron waves injected into the non-magnetic metal for different directions of the spin 

Here, d is the angle between wave vector k^. and the normal to the surface 2 = 0, pointing 
into the sample. The total effective polarization P^^ of the current depends on the difference 
between the probabilities of tunnelling for different cr. 

The functions \E''^'')(ri; ai) for each spin direction have the same spatial dependence as those 
for a contact between non-magnetic metals. In following calculations we use the asymptotic 
expression for \E'*^°)(ri; o"/) valid for A;a <^ 1 [9| 

^^'\r;a)=tM'\r)^^, (6) 

where 

= %&e*' (l + -L], (7) 



2ikr'^ \ ikr 

ijj^ is the spinor 



^T=(Jl ^i=(?)^ (8) 



and k = \/2m*e/h is the magnitude of the electron wave vector k in the non- magnetic metal. 
Note that the wave function ^(°)(r) ([ZD is zero in all points on the surface z = 0, except the 
point r = (at the contact) where it diverges. This divergence is the result of taking the 
limit a — > in the integral expressions for \E'^°^(r) [sl, [2o|. Yet, Eq. ([71) gives a finite value 
for the total charge current through the contact as obtained over a half-sphere of radius r 
with its center in the point r = for r — 0. 

The solution to Eq. ([!]) in linear approximation in the potential U, Eq. ([2]), can be written 
as _ _ 

§(r;(T)=t^ ^(°)(r)^;^ + (^(^5^± Jcosa) + Jsina^_^j ^(^)(r) ; (9) 

where the sign (±) corresponds to a =t, |. We have introduced the notation 

_ 2m*k ~ m*k ^ . , 

9 = ^9. J=^J/ieff (10) 

for the dimensionless constants of interaction for potential and magnetic scattering, respec- 
tively. The wave function scattered by the cluster is given by 

^^'\r) = dr'G+(r,r')D(|r'-ro|)vI/W(r'), (11) 

which undergoes specular reflection at the metal surface. Aiming for a solution for the wave 
function \E' ^.(r) in first approximation in the small parameter \t^\ ^ 1 we substitute the 
electron Green function G'(]"(r,r') of the non-magnetic half-space in Eq. iQ, 

Go+(r, r') = G+o(|r - r'|) - G^,{\r -r'\), (12) 
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where r = (p, —z) is the mirror image of the point r relative to the metal surface, and Gqq 
is the Green function for free electrons, 

„ikR 

GooiR) = -^; R=\r-T'\. (13) 

The wave function (jH]) enables calculation of the charge current density j and spin current 
density j'-^-*, and the expectation value of the spin s. They are obtained as the sums of 
independent contributions of the two electron groups (/ = 1,2) 



jYl* /L^ \ ) ri=r2=r 



Z=l,2 



j'''W-^E[(v5i-.5-5^-.H]_; (15) 

1=1,2 

s(r) = V f^+CT^) . (16) 

The Eqs . f lT^ - f|T6|l with wave function ([9]) describe the so-called tunnelling contributions 
(see, 16[). They are proportional to the tunnelling probability and define the measurable 
quantities which can be obtained after averaging over the energies of the tunnelling electrons 
and wave vector directions (see. Sec. III). 

In absence of the cluster 8^ = Sy = 0, and the local magnetization So(r) due to itinerant 
spin polarized electrons is oriented along the 2;-axis. The spin polarization in zeroth approx- 
imation, S2o(r), which is calculated from the wave function ([6]), is a monotonic function of 
coordinates 

»=»w-(w=-fen(^)'(i+^). (17) 

The electron scattering by the spin-depended potential ([2]) changes the value and the 
direction of the vector So(r). Components and Sy appear due to electron scattering by 
the cluster and they are subject to quantum interference between transmitted and scattered 
waves. As a result of the interference the spin components perpendicular to the z-axis are 
oscillatory functions of the coordinates, while s^o acquires a small oscillatory component 
proportional to . 

{ ^(r) } = (I^tI' ± l^il') Jsina |vl/(0)(r)| |v^(^)(r)| | } (^,(r) - ^^W) (18) 

Here |\l/^*^(r)| and v'j(r) are the absolute values and the phases of the coordinate part of the 
wave functions ([9]) (see, Eqs. ([7]) and f|TT]) ). 

Figure 2 shows the spacial distribution of the a;-component of the normalized spin 
density Sa;(r)/co in the vicinity of the contact. The normalization constant Cq = 
(|t||^ -|- ItJ^) J (/ca)^ sina/lGvr. The figure demonstrates that the spin component ^^^(r) 
changes sign in the space of the normal metal. The sign of depends on the difference of 
the phases v^j(r). In Fig. 3 we present a vector plot of the x-component of the spin current 
density j*-^-*. An intricate image of the distribution in orientation of j^"^'-* is visible. Note the 
lines at which the direction of the vector j*^^^ is inverted. 
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X {1/k units) 



FIG. 2: Gray-scale plot of the spacial distribution of the x-component of the spin density, Sx{t)/co. 
The coordinate plane xz in real space crosses the contact and the cluster of the radius = k^^, 
with its center in the point tq = (0,0, 15) k~^. 



III. CONDUCTANCE OF THE CONTACT, SPIN CURRENT, AND SPIN- 
TRANSFER TORQUE 

The total current through the contact can be evaluated by integration of the charge 
current density j(r), Eq. (JT^ . over a half-sphere centered at the point contact r = and 
covering the contact at z > 0, and integration over all directions of the electron wave vector 
on the Fermi surface e = ep . In the Ohm's law approximation and at zero temperature the 
total current through the contact / can be written as 

I = eVp{eFy j ^1 dQe{z)Q{k,)jr{r), (19) 

£=ep 

where dQ and dfik are elements of solid angle in real and momentum space, respectively, 
pi^ep) is the electron density of states at the Fermi energy, ep, for one spin direction, fc^ is 
2;-component of the wave vector, j.r(i") is the radial component of j(r) iHM . and Q{x) is the 
Heaviside unit step function. After integration of Eq. f[T^ the conductance G of the contact 
takes the form 

1 + I (g + Pes J cos aj W{ro) 



G — — — Gq 



(20) 

e=e-p 

where Go is the conductance of the contact in absence of the cluster 

Go = {4, + 4,] ff}':"L (21) 

727r {m*Uo) 

Here, fepo- and kp are the absolute values of electron wave vectors at the Fermi level in 
magnetic and non-magnetic metals, respectively; Peff is the effective spin polarization of 
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FIG. 3: Vector plot of the direction of the j;-component of spin current density, j*-^-*. The contour 

ne : 

chosen ro = and tq = (0,0, 15) k^^ . 



hnes correspond to = 0. The plane xz crosses the contact and the cluster; as in Fig. 2 we have 



the current injected through the contact (see Eq. ([5])); the constants g and J are given by 
Eqs. ([TO]), and 

W^(ro) = j dr'D(lr'-rol) yi{kr')n{kr'), (22) 

where ji{x) and yi{x) are the spherical Bessel functions. Equation (!20|) coincides with the 
result for a tunnel point contact between non- magnetic metals ^ when Peff = and fcr^ -C 1. 
When the radius of action td of the function D(|r — ro|) is much smaller than the distance 
between the contact and the cluster center ro, W{yq) is an oscillatory function of kr^ for 
krD > 1, but the oscillation amplitude is reduced as a result of superposition of waves 
scattered by different points of the cluster. The integral W^tq), (12^ can be calculated 
asymptotically for tq 3> ro, kvQ ^ 1, and fcrn ^ 1. For a homogeneous spherical potential 
L>(|r|) = V^^6(rD — r) {Vd is the cluster volume) the function W^Tq) takes the form 

iy(ro)-^f^V^^^-#, (23) 

where d = 2rD is the cluster diameter. The last factor in Eq. ( [23i) describes the quantum 
size effect related with electron reflections by the cluster's boundary. Such oscillations may 
exist, if the cluster boundary is sharp on the scale of the electron wave length. Fig. H] shows 
the dependence of the amplitude of the conductance oscillations on the cluster diameter. It 



8 




FIG. 4: Dependence of the the oscillatory part of the conductance on the tip position on the metal 
surface for different cluster diameters. The pg-coordinate is measured from the point Pq = at 
which the contact is situated directly above the cluster; tq = (0,0,10)/^^?; g = 0.5; J = 2.5; 
Peff = 0.4; a = 0. 



demonstrates that a 7r-phase shift may occur resulting from interference of electron waves 
over a distance of the cluster diameter. 

In Eq. f l20|) the term proportional to Pcfr takes into account the difference in the prob- 
abilities of scattering of electrons with different a by the localized magnetic moment /igg. 
It depends on the angle a between the tip magnetization and ^tgfj as cos a. The same 
dependence was first predicted for a tunnel junction between ferro magnets for which the 
magnetization vectors are misaligned by an angle a , and this was observed in SP-STM 



experiments [10 



Similar to the electrical conductance 
can be calculated 



20|) . the total spin current for each spin component 



GoV 



PeS + ^{ Pcsg + J COS a] W{ro) 



-I e=ep 



QGqV 



sm a 



evr 



(24) 



(25) 



J £ = £-p 



For our choice of the vector /Xgg, J'-^-* = 0. The value of the z-component of the spin current 
/(^) ([23]) is determined to a large extent by the polarization Pefr dS]) of the initial current. The 
oscillatory part of I^^^ is modified by the addition of a term due to spin-flip processes on the 
cluster. The spin current component perpendicular to initial direction of polarization, J^^\ 
has only a term that oscillates with rg, and which disappears when the magnetic moment 
^gfj is aligned with the z direction. 

The spin-transfer torque T acting on the magnetic moment /Xgfr is given by 



J 

h 



rfr'D(|r'-ro|)/Xeff X (s(r')); 



(26) 
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where 



(s(r)) = eVp{ep) 



-Sir 



(27) 



e=ep 



s(r) is defined by Eq. ( fT6ll . This torque is related with the spin polarized electron tunnel 
current. In linear approximation in J only the spin density contribution Szo{i^) f|T7|) . without 
interaction with the cluster, should be taken for the calculation of the torque fl26|) . In this 
approach = = 0. For large tq ^ we obtain 



T 



3GnV 



sm a- 



evr 



PpffJ 



1 + 



{kroY 



(28) 



e=ep 



The dependence of Ty ( I28l) on the angle a agrees with the dependence of the spin transfer 
torque in tunnel junctions between two ferromagnets with different directions of magnetiza- 
tion (l6|,[2l|. 

In this paper we do not aim to investigate the dynamics of the cluster magnetic moment. 
We only note that once the spin-polarized current-induced torque pulls the magnetic moment 
away from alignment with Ha, the cluster moment will start precessing around the field axis. 
The Larmor frequency is defined by the magnetic field due to combining the external field 
Ha and the effective field produced by the polarized current H^ff — —J (s(rQ)) /(7c/^b {Qc is 
the cluster 'g-factor'). The precession of the cluster magnetic moment gives rise to a time 

The study 



modulation of the SP-STM current as for clusters on a sample surface [22|, [23 



of non-stationary effects provides a further means of obtaining information on the cluster 
and the spin polarization of the current inside the sample. 



IV. DISCUSSION 

In summary, we have studied theoretically the current and spin flows through a tunnel 
point contact between magnetic and non-magnetic metals when the tunnel current is spin 
polarized in the geometry of SP-STM, Fig. 1. Electron spin flip processes due to a magnetic 
cluster situated in the non-magnetic metal have been taken into account. These processes 
result in the appearance of components of the spin current density j^*-'(r) perpendicular 
to the source direction (taken along the contact axis z), and a finite expectation value of 
the spin s(r). We have analyzed the contribution of tunnelling electrons to the spacial 
distribution of j(^''2^)(r) and Sx,y{r). It is found that these are non-monotonic functions of the 
coordinates and undergo strong spacial oscillations originating from quantum interference 
between partial waves transmitted through the contact and those scattered by the cluster 
(see. Figs. 2, and 3). Specifically, between the contact and the cluster there are neighboring 
regions in which j^^^(r) flows in opposite directions (Fig. 3). 

The oscillatory correction, AG, to the conductance Go of the ballistic tunnel point contact 
strongly depends on the magnetic scattering ( [201) . 

AG ~ (g + PesJcosa^ sm2kro, kro > 1- (29) 

The effective spin polarization Pgfr of the tunnelling electrons, and the dimensionless con- 
stants of potential scattering g and exchange scattering J are given by Eqs. ([5]) and ffTOl) . 
Generally, for a single magnetic defect, which has a magnetic moment of the order of one 
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Bohr magneton, /ig, the spin part of electron scattering gives only a small contribution to 
the electrical resistance. For a magnetic cluster with /ig^ ^ /ig, the exchange energy can be 
larger than the energy of spin-independent interaction {J > g, see Eq. fllOp ). An interesting 
phenomenon may be found when PesJ > 5'- A change of the direction of the vector /x^g 
(the angle a) controlled by an external magnetic field is predicted to lead to a change in the 
amplitude of the oscillations, and for certain directions the amplitude may vanish, Gosc = 0. 
This large magneto-orientational effect provides a new mechanism for obtaining informa- 
tion on magnetic particles buried below a metal surface. Note that this effect is observable 
only for a spin-polarized current: if = t|, in linear approximation in J the changes in 
the scattering amplitude for spin-up and spin-down electrons balance each other and the 
magneto-orientational effect is absent. 

As a consequence of spin flips due to the interaction of the electrons with the cluster 
the oscillatory part AI^^^ of the spin current in the original z direction obtains a correction 
which depends on the exchange constant J and the orientation of the magnetic moment (12^ 

AJ(^^ ~ (^Pcsg + J cos sin2fcFro, kpro > 1. (30) 

A component of the spin current perpendicular to the z direction, \ is formed only by 
scattered waves and as the result of quantum interference it becomes an oscillatory function 
of the distance ro fl25l) 

/^''■^ ~ Jsinasin2A;F''^o, fcp'^o > 1- (31) 

The spin currents (130|) and (13T|) appear even in the case of non-polarized current through 
the contact and they are related to the magnetic scattering by the cluster. 

The torque, which acts on the magnetic cluster due to spin-polarization of electric current 
is a monotonic function of the distance ro (128!) to linear approximation in the exchange 
constant _ 

Tj^ ~ sinaPeff-'' (^o/^F^^o) ' ^f'^o > 1, (32) 

where zq is the depth of the cluster below metal surface. 

These results may be exploited in future experiments for detecting and investigating indi- 
vidual magnetic clusters buried below the surface of a host metal. Specifically, a comparison 
of the amplitude of the conductance oscillations for different directions of the cluster mag- 
netic moment allows determination of the exchange energy J/^efr and for a known value of 
the exchange integral J to find /Xg^. 
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